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Abstract 

Using a data sample of 2.59 x 10 7 ip(2S) decays obtained with the CLEO-c detector, we perform 
amplitude analyses of the complementary decay chains tp(2S) — > 7Xci! Xci ??7r + 7r~ and ip(2S) — > 
TXci; Xci — ► ri'n + n~ . We find evidence for an exotic P-wave rj'ir amplitude, which, if interpreted as 
a resonance, would have parameters consistent with the 7Ti(1600) state reported in other production 
mechanisms. We also make the first observation of the decay ao(980) — > r/'n and measure the ratio 
of branching fractions S(a (980) ->• rfir) /B(a (980) ->• r/vr) = 0.064 ± 0.014 ± 0.014. The tttt 
spectrum produced with a recoiling 77 is compared to that with 7/ recoil. 



* Present address: Pacific Northwest National Laboratory, Richland, WA 99352 
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I. INTRODUCTION 



Hadronic charmonium decays, in which charm and anti-charm quarks annihilate into 
gluons, provide an excellent opportunity to study light mesons. The combination of well- 
defined initial states and the availability of a wide variety of final states allows for the 
strategic selection of reactions to isolate and study different light meson systems. The 
decays Xci t]^ + ^ and Xci — ► ri'n + ii~ , in particular, have two interesting characteristics. 

First, since the quark and SU(3) flavor of the rj and rf are relatively well-known, one 
could, in principle, model these decays using what is known about the OZI rule and SU(3) 
flavor symmetry and learn about the tttt isoscalar states recoiling against them. Such a 
technique has been proposed for other Xcj decay channels []]. 

Second, the decays Xci — > riit + it~ and Xci — >• r(ix+ix~ provide an opportunity to search 
for exotic J PC = 1 h states in the t]tt and rj'n systems. In fact, the only two-body S"-wave 
decays of the Xci available in these channels necessarily have the rjir or rj'n system in a 
configuration with J PC = Two exotic candidates, the 7Ti(1400) and the 7Ti(1600), have 
been reported in other production mechanisms to have decays to t]tt and 7/71", respectively. 
The 7i"i(1400) has been reported primarily decaying to r]ir j2H5], while the 7Ti(1600) has been 
reported to decay to 77V [HI E], &i7r [SJ |H], /i7r [TO], and pn [TT] . 

We present amplitude analyses of the processes ip{2S) — > "/Xci] Xa — > ri-n + n~ and 
ip(2S) — > JXci] Xci Vj-n^-n^ in which we study the various tjtt and nir resonances pro- 
duced in the decays of the Xci- In Section [TT], we describe the data-selection process that 



results in a sample of Xci decays with estimated backgrounds below 5%. In Section III 
we describe our construction of amplitudes using the helicity formalism. Here we assume 
that the Xci i]7i + 7i^ and Xa — > ri'ii + 7i~ decays proceed through a sequence of two-body 
decays, where the intermediate states have well-defined quantum numbers. We pay special 
attention to the treatment of the nn S-wave, which utilizes independent experimental data 
on Tin S-wave scattering [12]. We also describe our fitting procedure, based on the unbinned 
extended maximum-likelihood method. 



The highlights of our analysis, detailed in Section IV 



are 



• evidence for a P-wave 77V amplitude, which, when parametrized by an exotic J PC = 
1 h resonance, has properties consistent with those of the 7Ti(1600) reported in other 
production mechanisms; 

• the first direct observation of the decay cto(980) — > rfir, a measurement of the ratio of 
branching fractions £>(a (980) — > 7/71-) //3(a (980) — > r/n), and a characterization of the 
ao(980) lineshape; and 

• the observation of qualitative differences in the hit system when it is produced against 
the rj or rj'. 

Finally, in Section |Vl we evaluate and discuss systematic errors. 



II. DATA SELECTION 



We select candidate events of the form ip(2S) — > jXa] Xa — > rj-n + Ti~ and ip(2S) — > 
TXci! Xa — > ri'-n + it~ using 25.9 x 10 6 ip(2S) decays collected by the CLEO-c detector at 
the Cornell Electron Storage Ring. We reconstruct the rj (rj ; ) in three (six) different decay 
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TABLE I. The decay modes of the 77 and rf that are used to reconstruct ip(2S) — > 7Xci, Xd 
7/7r + 7r~ and tp(2S) — > JXcl'i Xd — > rj'iT + -K~ , the branching fractions B of each decay mode [14J, and 
the final state topology reconstructed with the detector. 



rP Decay Mode B [%] Final State 

77 -)• 77 
7/ —7- 7T + 7r _ 7r 
r? — > 7r 7T 7r 



39.3±0.2 


37l(vr- 


-7T 


22.7±0.3 


37 2(vr- 


-7T 


32.6±0.2 


77l(vr- 


-7T 


17.1+0.3 


37 2(vr- 


-7T 


9.9±0.2 


37 3(vT 


-7T 


14.1±0.2 


77 2(tt- 


-7T 


29.3±0.6 


272(71-" 


-7T 


8.5±0.3 


771(tt- 


-7T 


4.9±0.2 


77 2(vr- 


-7T 



77 — > TT^TT~rj; 77—7-77 

77' — > 7r + 7T — 77; 7/ —7- -/T + 7T~7r 

77' —7- 7r + 7r — 77; 7/ — > 7r 7r 7r° 

77' — > 77T + 7T _ 

77' — > 7r°7r°77; 77 — > 77 

77' —7- 7r°7r°77; r/ — )■ 7r + 7r _ 7T° 



topologies comprising 94.6 ± 0.7% (83.8 ± 1.8%) of its total decay width (Table |T|). We then 
select the Xd using the energy of the photon from ip(2S) — > 7Xci- 

Final state photons and charged pions are measured by the CLEO-c detector [13] , which 
covers a solid angle of 93%. The detector has a 1 Tesla superconducting magnet enclosing 
two drift chambers and a ring imaging Cherenkov (RICH) system for tracking charged 
particles and particle identification. Enclosed inside the solenoid are also a barrel and two 
endcap Csl-crystal calorimeters. The energy resolution for 100 MeV (1 GeV) photons is 
5.0% (2.2%), while the momentum resolution for charged tracks in the drift chambers is 
0.6% at 1 GeV/c. 

Charged tracks are required to have momentum p > 18.4 MeV/c and originate within 
a cylindrical volume, with 20 cm length and 2 cm radius, centered around the interaction 
point. The ^ candidates are then required to have ionization losses (dE/dx) within 3cr of 
those expected for charged pions. Photons reconstructed inside the calorimeters, with polar 
angles | cos^l < 0.79 and 0.85 < | cos^l < 0.93, are required to have energy E > 20 MeV and 
separation from charged tracks. Two-photon four-momenta are kinematically constrained 
to select the 7r° — > 77 and 77 — )• 77 candidates, with a requirement that the respective 
invariant mass is within 3a of the 7r° or 77 rest mass. Finally, the total four-momentum of 
all of the final state particles of a given topology is kinematically constrained to the initial 
ip(2S) four-momentum and the x 2 °f the resulting fit is required to satisfy x 2 /d-o.f. < 5. 
If multiple combinations of tracks and showers within an event pass all of these selection 
requirements (which occurs for 1.6% of all selected events), only the combination with the 
best x 2 /d-o.f. is accepted. 

To select the decays 77 — > 7r + 7r~7r° and 77 — > 7r 7r 7r°, the invariant mass of the three pions 
must be between 540 and 555 MeV/c 2 (see Fig. [TJ. Similarly, to select the 77' in its various 
decay modes the invariant mass of its decay products must fall between 950 and 965 MeV/c 2 
(see Fig. [2]). If multiple combinations of particles within an event can be used to form an 
77 (occurring in 0.9% of selected events), one of these combinations is chosen randomly. If 
there are multiple 77' candidates (occurring in 1.4% of selected events), the event is discarded. 
Since the decays 77' — > ix+ix~r}\ 77 — > 7r 7r 7r° and 77' — > 7r 7r°77; 77 — > 7r + 7r~7r share the same 
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FIG. 1. The invariant mass of rj candidates after selecting a Xci candidate and applying background 
suppression criteria. Candidates in (a) are selected by requiring individual photon pairs be no more 
than 3(7 from the nominal rj mass; the arrows in (b) and (c) indicate the region used to select the 
rj candidates. 



final state topology, no requirement is made on the internal r\ mass. 

Specific backgrounds are further suppressed based on studies using a Monte Carlo (MC) 
sample of inclusive if)(2S) decays. For the '0(25') —¥ "/Xci] Xa — > r\ix+ix~ decay chain, the 
dominant backgrounds are due to ip(2S) —> r/J/ip and ip(2S) — > jXci] Xa — > iJ/^, where 
J/ip — >■ /x 4 ~yT ' . The first of these is suppressed by requiring the mass recoiling against the 77 
be separated from the J/ip mass by at least 20 MeV/c 2 . The second is only a background 
for the 77—7-77 mode and is similarly suppressed using the masses recoiling against the 77 
combinations, which are required to be more than 35 MeV/c 2 from the J/ip mass. 

The largest backgrounds in the ij){2S) —> JXci] Xci — > rfn + iT~ decay chain occur in 
the rj — ¥ 77r + 7r~ mode. We suppress J/ip backgrounds, as above, by requiring the masses 
recoiling against the 77 and 7r + 7r~ systems to be more than 20 MeV/c 2 away from the J/ip 
mass. In addition, we treat 7r + 7r~ combinations as and require their invariant masses 

be more than 15 MeV/c 2 away from the J/ip mass. There is also a substantial background 
from ip(2S) —¥ 7r°2(7r + 7r~), which we reduce by requiring that no two showers in an event 
be consistent with the 7r° mass within 3a. We also enhance the signal to background for the 
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FIG. 2. The invariant mass of rf candidates after selecting a Xci candidate and applying background 
suppression criteria. The arrows indicate the region used to select the rf candidates. 



rf — > 77T + 7r mode by making a loose requirement that the 7r + 7r invariant mass be between 
335 and 895 MeV/c 2 , which is motivated by the apparent p dominance in the n + ir~ system. 

One additional background for the ■^(25') — > -fXa', Xa — > ri'-n + it~ decay chain with 
rf — > 77r + 7r~ is from ip(2S) —> 7Xco; Xco 2(7r + 7r~) where the radiated photon converts 
to an e + e~ pair outside the tracking region. This is suppressed by requiring that the total 
energy of the two resulting showers is not consistent with the energy of the photon from 
ip(2S) — > 7Xc0) 2.e. ; not between 225 and 295 MeV, and that the cosine of the angle between 
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FIG. 3. The invariant mass distributions of the (a) r]Tr + ir~ and (b) rfn^ir - candidates from selected 
ip(2S) — > jt]tt + tt^ and ip(2S) — > r y7]'-K + iT~ decays, respectively, after all background suppression 
criteria have been applied. The solid arrows indicate the regions used to select the Xci signals. 



the two showers is less than 0.97. 

Figure [3] shows the invariant mass distributions of (a) the r]iT + 7T~ and (b) the r]'-K + iT~ 
candidates after combining all of the decay modes of the 77 and 77'. We select the Xci 
by requiring that the energy of the photon radiated from the ip{2S) be between 155 and 
185 MeV (indicated by the arrows in Fig. [3]). Our final data samples consist of 2498 and 
698 events in the ip(2S) —> JX&] Xci i]ii + it~ and i/j(2S) —> JXd] Xci ri'-n + -n~ decay 
chains, respectively. The background is estimated by fitting the data in Fig. [3] using a reverse 
Crystal Ball shape [12] to describe the signal. The background and Xc2 peak are described by 
a second order polynomial and a double Gaussian, respectively. Peaking backgrounds have 
been subtracted by fitting the Xc candidate mass distribution in 77 W mass sidebands; such 
backgrounds are negligible in all cases except the 77' — > 77r + 7r~ decay mode. The estimated 
signal purity for the 777r + 7r~ (r/V+Tr - ) decay channel is 97.5% (94.6%) with an uncertainty 
of 0.3% (1.3%). 



III. AMPLITUDE ANALYSIS 



We perform amplitude analyses to disentangle the substructure present in the x& ~ > 
t]7t + 7t~ and Xci ~^ Vj-n^-n^ decays. We assume that the three-hadron decays of the Xci 
proceed through a sequence of two-body decays, where one participant is the "isobar," a 
bound state of either i]^^ or with total angular momentum J, and the other is a 

stable, non-interacting meson (the 7r T or 77 W) produced with an orbital angular momentum 
L with respect to the isobar. All possible x& decays through isobars with J < 4 are listed 
in Table M 

The general idea of an amplitude analysis is to fit the distribution of events observed 
with the detector to a coherent sum of physically-motivated amplitudes that describes the 
dynamics of the intermediate states. We can define /(x), the number of observed events per 
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TABLE II. A list of Xci decay modes for all possible isobars with J < 4. 



Xci Decay Mode 


L 


Isobar J 


a vr; a -> 




P 


0++ 


7T17T; 7Tx — 


• j/'V 


S,D 


1-+ 


027T; 02 — > 




P,F 


2++ 


047T; 04 — > 




F,H 


4++ 


/or/ (,) ; /o 


— )• 7T7r 


P 


0++ 




—> 7T7T 


P,F 


2++ 




—> 7T7T 


F,H 


4++ 



unit phase space, as 



IU) 



M^,,A 7 



E 



^ / M v ,,A 7 ^M v ,,A 7 ( X y 



(1) 



where a indexes the x& decay amplitudes and and A 7 index the polarization of the 
ip(2S) and the helicity of the photon, respectively. We use x to denote a set of kinematic 
variables, e.g., angles and invariant masses, that provide a complete description of the event. 
The value of the decay amplitude at a location x in this multi-dimensional space is given 
by Aff a 7 ( x )- The real fit parameters Vfo A ^ determine the relative strengths of each Xci 
decay amplitude. 



Section HI A discusses the construction of the decay amplitudes used in the fit. Sec- 



tion III B discusses the application of the extended maximum likelihood technique to this 



analysis in order to determine the optimal values of A that describe the data. 



A. Amplitude construction 



1. General amplitude structure 

The amplitude for a given Xci decay mode a depends on the set of observed final state 
event kinematics x, the assumed polarization of the initial state ip(2S), denoted M^, and 
the helicity of the final state photon A 7 . The general form is given by 

A X =±1,0 M' X =±1,0 

(1M' X \LM' L , JM' J )Y^\Q' I , ^¥^(91 <f>' h )p L q J T a (s), (2) 

M' L ,M'j 

where summations in the second line are performed over all possible values M' L and Mj, the 
projections of L and J, respectively. We briefly provide a term-by-term description of this 
expression. 

The first factor in Eq. C(M^p, X 7 ,X X ), is used to transform the helicity amplitude 
for the radiative decay into the multipole basis. The ip(2S) —> 7Xci radiative transition 
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is dominated by the electric dipole (El) transition, while the magnetic quadrupole (M2) 
transition contributes ~ 3% [16] of the total rate. In our analysis, we use the El amplitude 
to derive our results and check the sensitivity of the results to the presence of a small M2 
amplitude. The El or M2 amplitude can be constructed with the following choice of C: 

C(M i> ,X J ,X x ) = y|^ W7 _^0 7 ,£ 7 ,O) x 

(£a 7 ,A x ,i - S Xl -i6x x ,-i + (5a 7 ,i - <5a 7 ,-i)<5a x ,o) for £71, or 
(<*a 7 Ax.i - <*a 7i -i5a x ,-i - (<*A 7 ,i - 5a 7 ,-i)^a x ,o) for Ml. 1 ' 
In order to describe the angular distribution of the final state particles we measure angles 
in two coordinate systems which are depicted in Fig. [5] and related in Eq. §2§ by the D- 
function at the end of the first line. The angles # 7 and 7 are the polar and azimuthal angle 
of the radiated photon in the ip(2S) rest frame, where z is given by the e + beam direction 
and y is (arbitrarily) defined as upward in the laboratory. (The amplitude is uniform in 7 .) 

The two spherical harmonics in the second line of Eq. ^ provide a description of the 
angular distribution for the initial xa decay and the subsequent isobar decay for various 
values of isobar-hadron orbital angular momentum L and isobar angular momentum J. The 
angles Q\ and 0' 7 are the polar and azimuthal angles of the isobar in the x c i-helicity frame, 
defined as the rest frame of the Xci with z'-axis along the photon momentum and y'-axis 
perpendicular to the plane formed by the ip(2S) and photon three-momenta. Finally, the 
angles 9' h and 4> ; h are the polar and azimuthal angles of h, one of the hadrons produced in 
the isobar decay, after boosting the momentum of h in the Xci-helicity frame to the isobar 
rest frame [TTj . All values of M' L and Mj are summed with appropriate Clebsch-Gordan 
coefficients to create the initial xa state with one unit of total angular momentum and z' 
projection M' x 

The "breakup momentum" in a decay of 1 — > 2 particles is given by the momentum 
of one of the produced particles in the rest frame of the parent. We denote the breakup 
momentum of the initial Xci decay and the secondary isobar decay by p and q, respectively. 
Finally, the term T a (s), described in detail in the next section, is a function of the invariant 
mass squared of the isobar and describes the two-body dynamics in the decay. 

To impose isospin symmetry in the decays xa — > ajvr T , we write the ajn amplitude as 

1_ ( .at 

V2 



^mXa 7 ( X ) — ~/K ( ^A^,A 7 ( X ) + ^M V ,,A 7 ( X ) ) J (4) 



where the distinction between the two terms on the right-hand side is the interchange of the 
7r + and 7r~ four-momenta in the calculation of the relevant kinematic variables. A similar 
symmetrization is used in the construction of the ttitt amplitude. 



2. Two-body dynamics 

We use three different formulations of T a (s) [in Eq. §2§] to describe the isobar decay 
amplitude and phase as a function of s, the invariant mass squared of the isobar decay 
products. For all intermediate states except the a (980) and the nn S- wave we use a Breit- 
Wigner distribution, 

T *{s) = — 2 : TTTT' ( 5 ) 
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with 

2 



Bj(q(s)) 



(6) 



Po L Bj(q ) J 

where mo and To are the isobar mass and width, respectively. We define the breakup 
momentum q Q = q{m%). Likewise, the available phase space is given by p(s) = 2q{s)/y/s, 
and po = p( m o)- These factors are used in conjunction with Bj(q), a spin-dependent Blatt- 
Weisskopf barrier penetration factor [TH] , to construct the mass-dependent total decay width 
given in Eq. (|6|. 

To describe the ao(980) line-shape we use a three-channel Flatte formula [19J. In addition 
to the common decay modes a (980) — > r\ix and a (980) — > KK, we include a third decay 
mode, a (980) —> tj'tt, to provide a consistent description for both the Xci — > W^ 11 ^ an d 
Xd i]' 7l+7l ~ data. The parametrization takes the form 

ml- s-iJ2 c 9cPc 

where m is the ao(980) mass and g c represents a coupling to one of the ao(980) decay modes: 
rjTc, KK, or 77V. The factors p c are the phase space available for each of the three different 
final states. Following the technique suggested by Flatte to preserve analyticity at the KK 
and i]'tt thresholds, we allow the phase space factors to become imaginary when s is below 
threshold for a particular decay channel [T9] . 

For the inr S'-wave, we utilize an analysis of tttt scattering data [12] that provides two 
independent amplitudes for tttt — > tttt and KK — > tttt production mechanisms. Specifically, 
we attempt to model both direct production of Xci V 7171 with the tttt in an S'-wave and 
also the process xa t]KK — > t]tttt where the KK S'-wave intermediate state rescatters 
into tttt S-wave. These two amplitudes, labeled S n7T and Skk respectively, are constructed 
to be consistent with existing data in the region of tttt invariant mass below 2 GeV/c 2 , where 
the S-wave is expected to be significant. 

To account for the s-dependent differences between tttt scattering, from which the am- 
plitudes are derived, and production in Xci decay (see Fig. [5J, the S n7T scattering amplitude 
is rewritten in a form N(s)/D(s), and the numerator is replaced by the first two terms in a 
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(a) 




FIG. 5. A diagram of the mr scattering process (a) and the Xd decay process (b). Both processes 
may have intermediate pion or kaon loops. The dark grey interaction represents that obtained 
from scattering data, while the light grey interaction is some unknown s-dependent production 
amplitude in Xcl decay. 



series expansion 

SUs) = 1 -^f = Sl(s) + cSUs) (8) 
where the conformal transform z(s) is given by 

z(s) = Vj+j^v ffSEI (9) 
V s + s o + V - s 

The parameter so, used to adjust the left-hand cut in the complex plane, is set to so = 
1.5 GeV 2 /c 4 . Since the Skk amplitude only peaks in a narrow region of s, it is assumed 
that a similar s-dependent modification of the production amplitude is not necessary. The 
magnitudes and phases of the Skk, S® n , and amplitudes used in the fit are shown in 
Fig.© 

Given the definitions above, we can write an expression for the dynamical portion of the 
mi S-wave amplitude in the fit: 

T ( ^ )s (s) = Sl(s) + cSUs) + kS KK (s), (10) 
where c and k are real parameters that determine the relative sizes of the components. 



B. Fitting technique 

The extended maximum likelihood method is utilized to determine the best model pa- 
rameters to describe the data. The likelihood can be written as 

r- e ~^ N T\ C(x,)/(x.) 

N\ ll/C(x)J(x)dx' 1 j 
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FIG. 6. The (a) magnitudes and (b) phases of both terms of the irn — > irir (dashed and dotted, blue 
online) amplitude and the KK — > tttt (solid, red online) amplitude as a function of the two-pion 
invariant mass squared. 



where x, as above, is a position in the multi-dimensional space that spans the event kine- 
matics. The functions £ and / describe the efficiency of the analysis criteria and the model- 
predicted density of events in this space, respectively. The total number of observed events 
is N, and the model-predicted number of events is represented by 

p. = J C(x)/(x)dx. (12) 

In practice one varies the fit parameters to minimize the function 

-21n£ = 2^y C(x) I(x)t2x - In /(x*) j + K, (13) 

where k is a constant term that is not included in the minimization procedure. The integral 
on the right hand side of the equation is evaluated using MC techniques to determine the 
average value of C(x)/(x). Namely, we generate an MC sample of N g signal events that are 
distributed uniformly in phase space. These events are subjected to our selection criteria 
and yield a sample of N a accepted events. For the current study N a /N m 60. From this 
sample we can compute 

r v Na 

/ C(x) J(x)dx = V (C(x)/(x)> 7 ( x *)> ( 14 ) 

" i=l 

where V is the volume of phase space that spans the event kinematics. 

In the context of the model discussed in the previous section, /(x) is as given in Eq. (fl]). 
Finally, after substituting Eq. (14) into Eq. (13), making a change of variables x ^ — > 



x 7 /vV, and collecting constant terms into k', we can write the expression that is mini- 
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mized by the fitter: 



-21n£ - k' 



N„ 



N, 



E E 



1=1 M^,\ry 



N 



i=l \ M^,,A 7 



E/ ^ / M V ,,A 7 ^M^,A 7 ( X '«) 
E/ ^H/> • A 7 • A 7 ^ X 



(15) 



Note that k' is a constant that only depends on V and the number of events N; therefore, 
its actual value is not needed to construct a ratio of likelihoods (£) for two different models 
that describe the same set of data. 

In the fit, we constrain = V{*_-y = V* t x = Vf^ _ x for each decay amplitude. Recall 
that the coefficients C defined in Eq. ^ ensure the proper linear combination of ip(2S) and 
photon helicities is used to generate the desired El or M2 transition amplitude. 

In order to compare the results of our unbinned fit with the data for a given kinematic 
variable, we take the signal MC sample that passes our event selection criteria and weight 
each entry in a histogram of the given variable by the function /(x), defined above. To isolate 
contributions from various Xci decay modes, we can restrict the sum used to compute the 
weight /(x). 

The results of the fit can be cast as an acceptance-corrected "fit fraction" for each Xci 
decay amplitude a given by 



T — 

J OL 



/E 



M^,,A 7 



4 C 



Ai^,,A 7 A^^,,A 7 



G?X 



//(x)cfa 



(16) 



The integrals of the amplitudes over phase space are evaluated using MC techniques, similar 



to Eq. (14), except the sum is performed over the N g generated MC events. Note that due to 
interference between the decay amplitudes, x (x), the sum of the fit fractions, ^ a J", 
is not constrained to unity. 

We compute errors on the fit fractions by propagating the errors on the fitted values 
of \ through Eq. (16). We stress this point because such errors are not suitable for 



estimating the statisticaTsignificance of amplitudes with relatively small fit fractions. We 
obtain the statistical significance of amplitude a by computing the ratio of the likelihood 
of the null hypothesis (J-" a — > 0) to the likelihood of our baseline fit, which is derived by 
retaining known or possible ij^n and titt resonances that have a significance greater than 
one standard deviation. 



IV. RESULTS 



Figures u(a) and |7Fd) display Dalitz plots from the x& ~^ V n+7T ~ an d Xa — > ri'-n+n^ 
decays. We perform the amplitude analysis, discussed in the previous section, independently 
on both the x& t]7i + ti^ and x& t]'7t + 7i^ samples. In both cases our signal MC sample 
includes all of the various r\ and rj' decay modes that we reconstructed in the data populated 
according to their known branching fractions [TJ]. Figures [7^b) andj^c), and Figs. [7](e) 
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FIG. 7. Dalitz plots and corresponding projections from the Xcl 
(d-f) decays, overlaid with the baseline fits (solid lines). 
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and [7[f ) show projections of the Dalitz plot for the x& "^n^ix and Xci — > T]' , k + 'k data 
samples, respectively, along with the baseline fit for each sample. 

The projections of data and our baseline fits, separated into the various amplitudes, 
are shown in Fig. |8} For those isobars with J > we try multiple values of L (see Ta- 
ble[n|; only in the case of f2(1270)r]^ do we find a statistically significant contribution from 
higher orbital angular momenta. We report the sum of these P-wave and F-wave contri- 
butions as f 2(l270)r]^ ; the individual P and F wave fractions, as expected, have strongly 
anti-correlated statistical and systematic errors. We also separate the various (nix) s ri^ con- 



tributions given in Eq. (10). In the fits to the Xci t]'ii + ix~ sample we fix the txtx >S-wave 
parameter c in Eq. ( 10 ) to zero; allowing this parameter to float yields a value that is statisti- 
cally consistent with both zero and the value obtained in the higher-statistics x& ~^ i]7T + 7T~ 
fits. In all cases the parameters describing the masses and widths of the intermediate res- 
onances are fixed in our baseline fits to enhance the stability of the fit. Both the 7Ti(1600) 
and /4(2050) parameters are fixed to values that maximize the likelihood. We systematically 
explore uncertainties on the parametrization of the amplitudes as discussed in Section M 



A quantitative summary of the baseline fits appears in Table III From the fit one can 



compute the total acceptance-corrected event yield in either the r/7r + 7r or 7]tt + tt final 



states, which is the denominator of Eq. (16). If we denote this quantity N(r]7r + n ) or 



N(r]'TT + -K ), respectively, then we can compute the branching fractions for Xci to these final 
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TABLE III. Summary of results of the baseline fits. The first and second errors are statistical and 
systematic, respectively. The third error, where reported, is from the external value of B(ip(2S) —> 
7Xci)- Amplitudes that are preceded by an asterisk (*) are not part of the baseline fits but have 
been included to determine upper limits. The listed fit fractions and significances (N a ) are obtained 
when the amplitude is added to the baseline fits. 



Xci Decay Mode 


T [%] 


B(Xd -»■ r/Wvr+vr-) x T [10~ 3 ] 






- 


4.97 ±0.08 ±0.21 ±0.22 


- 


a (980)vr 


66.2 ±1.2 ±1.1 


3.29 ±0.09 ±0.14 ±0.15 


> 10 


a 2 (1320)7r 


9.8 ± 0.8 ± 1.0 


0.49 ±0.04 ±0.05 ±0.02 


9.7 


(7r + 7r~)s?7 


22.5 ±1.3 ±2.5 


1.12 ±0.06 ±0.13 ±0.05 


> 10 




12.1 ±1.7 ±5.6 


0.60 ±0.08 ±0.28 ±0.03 


> 10 




3.4 ± 0.9 ± 1.5 


0.17 ±0.05 ±0.07 ±0.01 


6.0 


SkkV 


3.1 ±0.6 ±0.4 


0.15 ±0.03 ±0.02 ±0.01 


9.4 


/ 2 (1270)ry 


7.4 ± 0.8 ± 0.6 


0.37 ± 0.04 ± 0.04 ± 0.02 


> 10 


h (2050)7? 


1.0 ±0.3 ±0.3 


0.05 ±0.01 ±0.02 ±0.00 


5.2 


*7Ti(1600)7T 




< 0.031 


0.7 






1.90 ±0.07 ±0.08 ±0.09 




a (980)vr 


11.0 ± 2.3 ± 1.8 


0.21 ±0.04 ±0.04 ±0.01 


8.4 


a 2 (1320)7r 


0.4 ±0.5 ±0.6 


< 0.031 


1.4 


(ir + ir-)si] 


21.6 ± 2.7 ± 1.2 


0.41 ±0.05 ±0.03 ±0.02 


10.2 


SLrf 


7.0 ±2.2 ±2.3 


0.13 ±0.04 ±0.04 ±0.01 


6.6 


S K kv' 


8.4 ± 1.5 ± 1.3 


0.16 ±0.03 ±0.02 ±0.01 


7.5 


h (1270)7?' 


27.0 ± 2.9 ± 1.7 


0.51 ±0.06 ±0.04 ±0.03 


> 10 


*U (2050)7?' 




< 0.010 


0.4 


7Ti(1600)7T 


15.1 ±2.7 ±3.2 


0.29 ±0.05 ±0.06 ±0.01 


7.2 



states, B(xd — ► ??7r + 7r ) and B(xd — > 71+71 )> as 



pN(ri { ' 



N H2S) B^(2S)^ 1X ci)EiW 



(17) 



where N^s) is the number of initial ip{2S), 2.59 x 10 7 , and we use B(^(2S) — > JXd) — 
(9.2 ± 0.4) x 10~ 2 [H]. The sum over rj and rj branching fractions encompasses all r\ and 7/ 
decay modes in our signal MC sample, indicated in Table [1} The value p is the purity of the 
data sample in the %& region in Fig. [3j which is obtained as discussed in Section [TT} 

In what follows we discuss the results of the fits to each of the samples in detail, high- 
lighting the key results obtained from each fit. For each x& decay mode, we also compute 
the product B(Xci — > //'• ) 7r + 7r~) x J 7 , which can be interpreted as the branching fraction 
for the Xci decay to the isobar and spectator multiplied by the branching fraction for the 
isobar to decay to the 7] ( -'' ) 7i ± or final state. Dividing products with common factors 



(discussed in Section IV C ) yields Xci and isobar branching ratios. 
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FIG. 8. Invariant mass projections from the analysis of the Xci r/ir + ir (a,b) and Xci 7? / vr + vr 
(c,d) decays. 



Xci ?77r+7r decays 



The ?77r =t: and 7r + 7r _ invariant mass projections and corresponding amplitude contribu- 
tions from the fit to the i/j(2S) — > 7Xci; Xci — >• iiit + it~ sample are shown in Fig. [8^a,b) . The 
dominant amplitude in this data set is the cto(980), which, consequently, must be adequately 
parametrized to obtain a satisfactory fit to the data. To determine the ao(980) parame- 
ters, we exclude the data with 7r + 7r~ invariant mass below 1.7 GeV/c 2 , which removes any 
correlation with the 7r + 7r~ S-wave amplitudes. The fit to this restricted data set includes 
the a (980), a 2 (1320) and / 4 (2050) amplitudes, and we allow all four a (980) parameters 
to float. The resulting ao(980) parameters are given in Table IV, where the first error is 



statistical and the second error is systematic, obtained by trying various combinations of 
tttt isobars to fit the region in tt + tt~ invariant mass around 2.0 GeV/c 2 , a peak attributed 
to the /4(2050) resonance in the baseline fit. The ao(980) Flatte distribution parameters, 
which are consistent with a previous determination by CLEO [20] , are subsequently fixed in 
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TABLE IV. The values of ao(980) parameters compared to the previous CLEO analysis [20]. The 
first error is statistical and the second error is systematic, as explained in the text. 



Parameter [GeV/c 2 ] Ref. [20] [GeV/c 2 

mo 0.998 ±0.006 ±0.015 1.002 ±0.018 
g m 0.60 ±0.02 ±0.03 0.64 ± 0.05 
g KK 0.56 ±0.06 ±0.09 0.52 ±0.15 
g^n 0.00 ±0.15 ±0.07 



the baseline fit to the full data sample. It is worth noting that the do (980) lineshape in tjtt 
is rather insensitive to the a (980) — > r(ix coupling g v > n . In fact, the fit prefers a coupling 
of zero, but with large uncertainty. Our analysis of Xci — >• r(ix+ix~ data, presented later, 
directly extracts information related to this coupling constant. 

The 7T7T S- wave is parametrized as described in Eq. (10) with the parameters c and k 



floating in the fit. In Table |HI| we list the contributions of the three individual components 
of the 7T7T S"-wave. In principle, the magnitude and phase of the total tttt S-wave can be 
constructed by using the entries in this table to normalize three components depicted in 
Fig. [6} 

In order to fit the it + tt~ invariant mass distribution around 2.0 GeV/c 2 , we tried various 
known Tin resonances with J = 0, 2, and 4 and masses ranging from 2.0 to 2.3 GeV/c 2 . The 
best fit is obtained with a single spin-four resonance that has parameters consistent with the 
/4(2050) state listed by the PDG [14J. The mass and width of the /4(2050), as determined 
by our fit, are m = 2.080±0.025±0.010 GeV/c 2 and T = 0.160±0.035±0.040 GeV/c 2 . The 
systematic errors are obtained by varying the ao(980) parameters within the errors listed in 
Table IVl 



B. Xci rfir+ir decays 



The set of amplitudes used to fit the xa t(tx+tx~ data is listed in Table III, and the 
corresponding fit projections are shown in Fig. |8](c,d). The dominant isobar in this Xci decay 
mode is the / 2 (1270), coupled with the rj in both P and F waves. Two key results emerge 
from the analysis of the rfir spectrum. First, a P-wave intensity is necessary to describe 
the rfir^ mass spectrum in the region of 1.7 GeV/c 2 . Second, we find that a (980) — > ti'tt 
decays populate the rj'ir^ region near threshold. We discuss both of these findings in the 
subsections below. 



1. Evidence of the P-wave rj'ir amplitude 

A fit to the data without the xa — > ^iir amplitude is shown in Fig. |9](a). It poorly 
describes the data in in the ri'ir invariant mass region near 1.7 GeV/c 2 and greater than 
2.3 GeV/c 2 (due to the contribution from the isosp in-conjugate channel). Our baseline fit 
accounts for these deficiencies by introducing a Xci — > ^i^ amplitude where the 7Ti resonance 
shape is described by Eq. do]) and the parameters for the mass and width are determined 
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by the fit. Such a resonance has exotic quantum numbers, J PC = 1 h , and cannot be a qq 
meson. 

The statistical significance of the Xci ^i 71 " amplitude can be evaluated by examining 
the ratio of maximum likelihoods for the two fits. If we define 

_ C{ni-n excluded) , . 

C(ji\T^ included) ' 

then, in the absence of any true tt-^ti amplitude, —2 In A will approach, in the limit of infinite 
statistics, a x 2 distribution with the number of degrees of freedom equal to the number 
of additional free parameters in the fit that includes the ttitx amplitude. As indicated in 
the first line of Table [V] the value of —2 In A for the baseline fit with and without the ttitt 
amplitude is 53, and three additional free parameters are used to describe the ttiti amplitude. 
This results in a very small probability (« 10~ n ) that data are a fluctuation of the model 
used in the fit to Fig. [91(a), which does not have a nin amplitude. 

In order to search for and quantify the significance of other non-exotic alternatives to 
the 71"! 7r amplitude that is used in our baseline fit, we considered several alternative fits 
to the data, enumerated 2-6 in Table [V] and described briefly here. Fit 2 replaces the 
exotic amplitude with the a (1450), a state that has a known tj'tt decay channel. Fit 3 
adds an a2(1700) amplitude, which might also decay into the rfw final state. Fit 4 adds an 
/o(1710) state in the im channel; heavier im states of various spins produce no significant 
improvements in fit quality. Finally, fits 5 and 6 test the J = 1 assignment of the rfit state 
by attempting to fit the data with a new ao or a 2 r/'n state whose mass and width are 
floating in the fit. The fit returns a value of the mass and width of 2.5 GeV/c 2 and 1.4 GeV 
(1.6 GeV/c 2 and 0.1 GeV) for the new a (a 2 ) state, respectively. In the third column of 
Table |V| we list the change in — 21n£ from fit 1, the baseline fit that does not include an 
exotic amplitude. None of the alternate fits produces a change in likelihood as significant 
as the fit that includes the tcitt amplitude. Furthermore, we can test the significance of the 
7Ti7r amplitude in the presence of these alternate models by including this amplitude in each 
of the alternate fits. The resulting values for —2 In A are listed in Table [Vj Including the 
7Ti7r amplitude introduces just one extra degree of freedom since the mass and width of the 
7Ti are fixed to the values obtained in the baseline fit. The tcitt amplitude is least significant 
in the context of fit 5, but even in this fit the significance of the amplitude is 4.7 (\/22) 
standard deviations. 

It has been noted that intensity in the P-wave rj'n scattering amplitude does not have a 
unique interpretation in terms of QCD and/or hadronic degrees of freedom [2T]. A complete 
analysis of the amplitude and phase of the P-wave is needed. Unfortunately, it is impossible 
to extract this information from the data in a model independent way due to the relatively 
low statistics and three-body nature of this analysis. Furthermore, non-resonant rj'n P- 
wave interactions are not well constrained, which makes it difficult to systematically test 
the significance of an exotic resonance in the presence of a non-resonant P-wave background. 
In an attempt to probe the significance of the Breit-Wigner phase motion, we replace the 
7Ti Breit-Wigner parametrization in our baseline fit with an amplitude whose magnitude 
matches that of a Breit-Wigner function but whose phase is constant (independent of s) 
and a free parameter in the fit. In this fit, fit 7 in Table [V} we also float the mass and 
width of the Breit-Wigner shape that describes the magnitude of the amplitude. We obtain 
a mass and width consistent with those obtained by introducing a resonant ttitt amplitude 
and the improvement in the fit is 45 units for four additional parameters, which is not as 
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(a) Fit 1 
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(d) Fit 4 

Amplitudes used in Fit 3 

f : (1710)-,' 




(b) Fit 2 



■ Amplitudes used in Fit 1 
a (1 450)tc 




(e) Fit 5 



-Amplitudes used in Fit 1 
■ a (2500)7c 




(c) Fit 3 

Amplitudes used in Fit 2 

a 2 (1 700)n 




(f) Fit 6 



- Amplitudes used in Fit 1 
■a,(1 570)7t 
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M(ri'7t) [GeV/c 2 ] 

FIG. 9. Projections of the t/tt invariant mass for fits to the Xci v(ix Jr ix~ data when the exotic 
7Ti(1600) amplitude is excluded. Fits are numbered according to Table |V} In all figures the sum of 
all amplitudes in the fit is indicated by the solid black histogram. 



dramatic as is obtained by including a ix\tx amplitude. The significance of the tx^tx amplitude 
in the presence of this non-resonant P-wave drops to 3 standard deviations. Our choice of 
parametrization for non-resonant interactions in this test is somewhat arbitrary - other 
choices may yield variations in the significance of the resonant txitx amplitude. While this 
is suggestive of significant resonant behavior in the P-wave, such a definitive conclusion is 
difficult to make without a complete understanding of non-resonant P-wave interactions. 
We can, however, state that there is clear evidence for P-wave t/tt interactions. 

Finally, we examine possible correlations of the P-wave intensity with the parametrization 
of the txtx S-wave. An examination of Fig. [8](d) suggests that the region of txtx invariant mass 
below 1.0 GeV/c 2 may not be well-described in the baseline fit. A natural concern is that 
the significant txitx signal might be correlated with the txtx S-wave amplitude in this region. 
To address this concern we perform two additional sets of fits. First, we exclude all data 
with txtx invariant mass below 1.0 GeV/c 2 , remove the txtx S'-wave amplitude from our fits, 
and repeat the exercise outlined above to compare various alternate solutions without a tx\tx 
amplitude. We find that including the tx\tx still produces the most significant improvement 
in the likelihood. Second, we devise an alternate parametrization of the txtx S'-wave that 
is not derived from txtx scattering data. This parametrization utilizes a complex pole to 
describe the broad txtx peak often called the a in addition to Breit-Wigner resonances for the 
/o(980) and /o(1370), with parameters set to those provided by PDG averages [Hj. Such 
a technique of adding many resonances has been previously used to describe txtx S-wave 



interactions [22]. The fit with this alternate txtx S-wave parametrization is shown in Fig. 10 
and the fit quality is improved in the region of txtx mass below 1 GeV/c 2 . We then repeat 
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M(ti'ti) [GeV/c 2 ] M{nn) [GeV/c 2 ] 

FIG. 10. Invariant mass projections of a fit to the v'tttt data that uses an alternate parametrization 
of the 7T7T S-wave. (See text for description). 



fits 1 through 6 listed in Table [V] with this alternate nir S- wave substituted into our baseline 
fit; the results are qualitatively the same as those derived from our baseline analysis. The 
corresponding value of —2 In A for fit 1 is 56 units, and the least significant signal for the 
7ri7r amplitude is in the context of fit 4, where the significance is 4.4 standard deviations. 

To summarize, the best fit to the Xci t/tttt data is obtained when a tx\tx amplitude is 
included, where the tci is described by a Breit-Wigner lineshape with a mass and width of 
1670 ±30 ±20 MeV/c 2 and 240 ±50 ±60 MeV/c 2 , respectively. The significance of a P-wave 
tj'tt amplitude is greater than 4 standard deviations under all attempted variations of the 
fit, some of which are rather extreme and assume the existence of new conventional cto and 
a2 states. While our baseline fit assumes that the 1~ + rj'it amplitude can be described by 
a Breit-Wigner resonance, we cannot exclude other non-resonant P-wave rfir interactions 
that may mimic a tt\ resonance. Therefore, we conclude that evidence exists for a P-wave 
ri'ir scattering amplitude, which, if parametrized by a single Breit-Wigner resonance, has a 
mass and width consistent with the 71"! (1600) reported in other production mechanisms. 

Motivated by reports of a 7r 1 (1400) observed in the rjii spectrum [2HS], we test the sig- 
nificance of an additional P-wave r)7r and rfn resonance that has the mass and width of 
the 7Ti(1400) reported in Ref. [14J. In neither the rjinr fit nor the rj-n-n fit is the signal for 
such a state robust under variations analogous to those in Table [V] which were used to test 
the significance of the 7Ti(1600) amplitude. Averaging over the two charge conjugate decay 
modes, we obtain the following 90% confidence level upper limits: B(xci — > 7Ti(1400) ± 7r =F ) x 
S(7ri(1400) ± -> r/7r ± ) < 0.08 x 10~ 3 and B(x<* -> 7^(1400)^) x i3( 7 ri(1400) ± rfif^) < 
0.02 x HT 3 . 



2. Observation of a$ (980) — > rj'ir decays 



Figure 11 shows a projection of the rfir invariant mass for a fit that includes all amplitudes 
in the baseline fit except the ao(980)7r amplitude. This fit does a poor job of describing the 
data near the tj'tt threshold. Our baseline fit corrects this by introducing an ao(980)7r 
amplitude, where the ao(980) lineshape utilizes the same Flatte parameters for the mass 
and couplings to rjir and KK channels obtained from the analysis of the r]iT + TT~ data. The 
significance of the a (980)7r amplitude, when compared with this alternate fit, is 8.4 standard 
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FIG. 11. Projection of the rj'ir invariant mass for a fit to the Xd ~^ V 'Tr + ir~ data that excludes 
the ao(980) amplitude. The value of — 21n£ for this fit, with one fewer free parameter, is 70 units 
larger than the baseline fit. 



deviations. 

Like the studies above for the tcitc amplitude, we try a variety of alternate fits that 
include the ao(1450) and a2(1700) resonances in the rj'n channel and the /o(1710) and 
/2(1750) resonances in the nir channel to describe the intensity at the rj'n mass threshold. 
The maximum change in —2 In £ observed for any alternate fit was 16 units, compared with 
70 units when the a (980)7r amplitude is included. We can repeat this study with the t\\ix 
amplitude excluded from the fit and similar results are obtained. Finally, we also study 
the stability of the ao(980) amplitude when the alternate resonance-based parametrization 
of the nn S- wave, discussed above, is used. We find that the overall fit fraction of a (980) 
decreases to about 50% of the value obtained in our baseline fit; however, the signal remains 
significant at the level of 6 standard deviations. 



C. Branching ratios 

Using the product branching fractions for the two final states, we can construct two 
different types of branching ratios. If we divide products that contain a common tttt isobar, 
then we can compare the production of this isobar in Xci decay against an rj and if. If we 
divide products that contain a common rj^it isobar, the result is measurement of the ratio 
of rjTT to rfir branching fractions for that isobar. We discuss each of these types of branching 
ratios below. 

Since the rj and rj' have a well known composition in both the quark and SU(3) flavor 
bases, it may be valuable to examine ratios of branching fractions for Xci decay to tttt 



isobars and an rj or rj' . Figure 12 shows a comparison of the tttt spectra for the fits when 
the recoil particle is an (a) rj and (b) rj' . The two components of tttt S-w&ve corresponding 
to KK — > tttt (Skk) and tttt — > tttt (SVtt) scattering have been highlighted. There are 
qualitative differences in the two tttt S'-waves. Most notably, in the rj recoil case, the Skk 
amplitude dominates the SVtt amplitude producing a peak near the /o(980). However, for 
the rj recoil case, the interference between Skk and S nn produces a dip near / (980), which 
is more clearly visible in the coherent sum of these amplitudes depicted in Fig. [8Tb). A 
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FIG. 12. Comparison of the irn spectra and the tttt S-wave contributions to the Xci i]tt + tt^ (a) 
and Xci rj'Tt + -K~ (b) decays. 



TABLE VI. Branching ratios of the Xci to various tttt states recoiling against rj and r/.The errors 
are statistical and systematic, respectively. Limits are set at the 90% confidence level. Correlations 
in errors on the individual fit fractions have been accounted for in the error on the ratios. 



Xci Branching Ratio Value 

(™)sr//(mr)sv 0.37 ± 0.05 ± 0.06 

Skkv'/Skkt] 1.03 ±0.28 ±0.17 

S nn ri'/S nw ri 0.11 ± 0.04 ± 0.04 

f 2 (1270)777/2 (1270)?? 1.39 ± 0.20 ± 0.11 
/ 4 (2050)r?7/ 4 (2050)r7 < 0.20 



compilation of branching ratios for all tttt isobars is listed in Table VI 



Table VII lists the branching ratios for r/ir and r/'ir decay channels of various isobars. The 
fact that the ao(980) is below tj'tt threshold undoubtedly contributes to the small tj'tt to rjir 
branching ratio. The small a2(1320)7r amplitude in our baseline fit yields an upper limit 
that is consistent with previous measurements: £>(a2(1320) — > 7/7r)/£>(a2(1320) — > r]n) = 
0.037 ±0.006 [Hj. Finally we can place a lower limit on B(tti — > 7]'ti)/B('Ki — > rjir) by adding 
a 7Ti7T amplitude to our fit to the rjTTiT data. This limit indicates that the rj'n partial width 
of the 7Ti is much larger than the r)7r partial width. 



V. SYSTEMATIC ERRORS 



In general, the systematic uncertainties on the results of this analysis can be classified into 
two broad categories: those that affect the measurement of the fit fractions (J 7 ) and those 
that affect the measurement of the total branching fraction for Xci — >■ ri ( -^n + n~ . Table VIII 



presents a summary of the systematic errors on these quantities. Correlations are considered 
when assigning systematic errors to branching ratios and products. 

The total branching fraction of Xci — > i]^tt + tc~ is affected by uncertainties in the Monte 
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TABLE VII. Branching ratios and limits for the ao(980), a2(1320), and 7Ti(1600) decays into the 
rf'n and ijtt decay modes. The errors are statistical and systematic, respectively, and the limits 
are set at the 90% confidence level. Correlations in errors on the individual fit fractions have been 
accounted for in the error on the ratios. 



Branching 


Ratio 


Value 


a (980) -» 
a 2 (1320) - 
tq(1600) - 


r/Tr/oo (980) -»• 
-> 77^/02(1320) 
-> r/7r/7n(1600) 


t/tt 0.064 ±0.014 ±0.014 
-»• r/vr < 0.065 
->■ r]TT > 9.1 



TABLE VIII. Fractional systematic errors in percent on the measurements of Xci ~~ 
branching fractions and amplitude contributions. See text for explanations. 



Decay Mode 


N lp(2S) € 7T± 




E1/M2 


Bkg. 


T a (s) 


Total Systematic 




2.0 0.8 


3.4 


0.7 


1.2 




4.3 


a (980)7T 






0.2 


0.4 


1.3 


4.4 


a 2 (1320)7r 






0.9 


7.1 


7.6 


11.2 


(-K + 7T-)sV 






0.5 


3.4 


10.7 


11.9 








2.1 


2.8 


45.7 


46.1 








2.6 


7.3 


42.2 


43.2 


SkkV 






1.3 


4.2 


11.6 


13.1 


/ 2 (1270)r? 






0.9 


5.9 


6.2 


9.6 


/ 4 (2050)r? 






1.0 


27.8 


13.4 


31.2 




2.0 1.2 


2.9 


0.6 


2.0 




4.3 


a (980)7T 






1.0 


9.4 


13.1 


16.7 








1.6 


4.6 


2.6 


7.0 


°1T1T I I 






1.3 


21.3 


24.1 


32.5 


S K kv' 






0.8 


8.5 


12.4 


15.6 


72(1270)77' 






0.5 


3.2 


5.4 


7.5 


vri(1600)7r 






2.7 


10.3 


18 


21.4 



Carlo model of the track and photon efficiency, e n7V - and e 7 . We assume a 0.3% systematic 
error for each track and 1% systematic error each photon. The total error is obtained by 
assigning a systematic error to each decay mode of the 77 or rj' and then constructing a 
weighted average of these individual errors where the weights are given by the product of 
the branching fraction and detection efficiency, i.e., a weight proportional to the number of 
observed events, for each mode. Assuming that the systematic error is not dependent on 
location of the event in the r]^n + n^ phase space, these errors cancel in the determination 
of the fit fractions. 

Background events have the potential to affect our analysis in two different ways. First, 
our computation of the purity p of each sample obtained from fits to the spectra shown in 
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Fig. [3] may be subject to systematic bias that would affect the measurement of B(x c i 
r]^'\ + 7r~). Second, because the background level is small and its angular distributions are 
difficult to characterize, our fit does not include a background amplitude. It is assumed 
that the background distributes itself among the various amplitudes thereby leaving the 
fit fractions unchanged. We test this assumption and also probe the stability of our total 
branching fraction measurement by relaxing event selection cuts and repeating the analysis. 
Specifically, we relax the x 2 /d.o.f. requirement; widen the invariant mass regions used to 
select T], r]', and Xci] introduce some J/ip background by reducing the effectiveness of the J/ij) 
suppression criteria; and enhance the probability that photons reconstructed in the event 
are actually decay products of 7r° by reducing the effectiveness of the ir° veto requirement. 
We take the largest deviation from our baseline analysis as the systematic error due to 
background. 

The construction of our amplitudes assumes the radiative transition ip{2S) —> ^Xci is 
purely electric dipole El. However, the contribution of magnetic quadrupole M2 amplitude 
has been measured by CLEO [16J to be (2.76 ± 0.76)%. This slightly affects the angular 
distribution of the radiated photon and the polarization of the x&- We quantify the uncer- 
tainty due to this assumption by repeating the analysis using the measured CLEO value and 
assigning the deviation from our baseline analysis as the systematic error. The dominant 
effect on the measurement of B(x& T]('\ + 7r~) is due to the change in detection efficiency 
of the radiated photon, while altering the polarization of the Xci affects the fitted values of 
J 7 . While these two effects are independent, their source is fully correlated, and we take this 
into account when obtaining product branching fractions. 

The choice of parametrization for the two-body dynamics T a (s) has the potential to 
systematically bias the results. While the error in the total branching fraction due to such 
variations is negligible, the individual fit fractions can be strongly affected by variations in 
T a (s). We vary, individually, the mass and width of the 0,2(1320) and ^(1270) within one 
standard deviation as tabulated by the PDG [14J and repeat the analysis. We also vary the 
parameters of the ao(980) according to the uncertainties listed in Table IV The mass and 
width parameters of the 7Ti(1600) and ^(2050), which are fixed in our baseline analysis to 
their best-fit values, are also varied by one standard deviation, and the analysis is repeated. 
Finally, we test the sensitivity of our results to the parametrization of the titt 5-wave by 
scaling the magnitude of this distribution by a value that ranges from unity to 1.2 or 0.8 
linearly with s. That is, we try various random linear changes in the shape at the 20% level. 
In addition, we vary the parameter so in Eq. (J9|. The largest deviation from our baseline 
fits within this set of variations is taken as the systematic uncertainty due to the amplitude 
parametrization. 

Finally, we must use as inputs the number of ip{2S) events in the data sample and 
branching fractions for the relevant rj, rj', and ip(2S) decays. The number of ip{2S) decays is 
known to a precision of 2%. The other branching fractions and their errors are taken from 
the PDG review [H] and listed as a separate, external systematic error for each measured 
quantitiy. 



VI. CONCLUSIONS 

We present an analysis of ip(2S) — > 7Xci — > 77/ / ' ) 7r + 7r~ decays in which we study the pro- 
duction of various rj^ir and titt intermediate states. Both channels exhibit a signal purity of 
at least 95% and the majority of rj and rf decay modes are utilized in the analysis. The con- 
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tributions from various quasi two-body decays are extracted utilizing an unbinned maximum 
likelihood fit that spans the phase space of relevant kinematic variables needed to describe 
the decay. Two-body interactions in our model are parametrized by both Breit-Wigner and 
Flatte distributions, and we utilize a parametrization of the tttt S-wave interactions that is 
based on scattering data. The results presented here supersede those previously presented 
by the CLEO Collaboration [20]. 

We find evidence for an exotic tj'tt P-wave scattering amplitude at the level of 4 standard 
deviations under a wide variety of model variations. If we parametrize this amplitude as a 
Breit-Wigner resonance we obtain a mass and width that is consistent with the 7Ti(1600) 
state reported in the literature. While the best description of the data is achieved with a 
resonant tt\tt amplitude, it is impossible to exclude other mechanisms that contribute to 
the rfn P-wave amplitude. In addition, the Xcj r(ix^ix~ data provide the first direct 
evidence for the decay of ao(980) — > tj'tt. We measure the ratio of branching fractions for 
the tj'tt and iitt decay channels of the ao(980). The xa t]TT + Tr^ data allow us to extract 
the parameters of the do (980) lineshape in the context of a three-channel Flatte distribution. 
Finally we compare the tttt system when it is produced against an r] to that produced against 
an 7/ '. Our model for the tttt S'-wave interactions suggests that production via tttt — > tttt S'- 
wave scattering is suppressed with respect to KK — > tttt scattering when the system recoils 
against the 7/. We also extract similar branching ratios for the other tttt resonances used in 
the fit, the / 2 (1270) and / 4 (2050). 
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